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1 Introduction 

The mechanical properties of biopolymcrs like DNA and actin are well described 
by a model of semiflexible polymers, characterized by a tangent-tangent cor- 
relation length that is considerably larger than the monomer size. A proper 
mathematical description is the worm-like chain (WLC): an elastic space curve 
with bending modulus A, that directly relates to this correlation length or ori- 
entational persistence length P;, = A/ksT. 

In the living cell biopolymers are usually in a crowded, strongly confined 
environment. Understanding the essential difference in statistics between the 
two extremes of the mechanical, strongly confined limit and the critical point of 
the self avoiding random walk limit, was an important development in polymer 
physics. In the first limit, the so called Odijk regime, the confining potential 
dominates on length scales shorter than the persistence length while in the lat- 
ter case, the deGennes regime, the excluded volume interaction dominates the 
thermal fluctuations on length scales of the channel diameter. The general pic- 
ture is the following [16]. Starting from strong confinement in a narrow channel: 
the polymer has due to its stiffness no other option than to follow the channel 
monotonously. The length scale over which the polymer is unperturbed by the 
confining potential, the so called deflection or collision length, is the dominat- 
ing length scale. Upon increasing the channel size, or decreasing the potential 
strength, hairpins will form inverting the direction of the backbone of the poly- 
mer. The average distance between hairpins, the so called global persistence 
length, decreases with increasing channel size and becomes the length scale that 
determines the span of the polymer until this global persistence length becomes 
of the order of the oricntational persistence length. At that point the a transi- 
tion regime sets in in which the volume interactions, needed in the traditional 
de Gennes regime, are non existent relative to the thermal energy This is 
a consequence of the slender shape of the polymer. Finally for channel sizes 
considerably larger than the oricntational persistence length, when the channel- 
length over which the volume interactions have to be added to account for at 
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least 1 ksT shrinks to the diameter of the channel, the de Gennes description 
becomes applicable. 

The reduction in degrees of freedom caused by the confinement, can be 
expressed as an entropic cost that is often a decisive ingredient to statistical 
physics calculations, for systems in thermal equilibrium, but also for the dy- 
namics of the system. Examples are the ncmatic transition of semiflcxible poly- 
mers, the ejection of DNA from viral capsids [12], translocation of DNA through 
nanopores [14], the elongation and dynamics of DNA in nanochannels [19] and 
the dynamics of actin networks [G] . From an application point of view the possi- 
bility to lay out DNA almost fully stretched in a nanochanncl makes it possible 
to study DNA and its interaction with proteins on the base pair level [18]. 

In many cases the WLC model is not sufficient though and it becomes crucial 
to understand how twist effects change this picture. This is for example the case 
in understanding the physics of transcription, the dynamics of supercoiling of 
DNA in vivo and the torque molecular motors can produce. In this paper the 
influence of torque on the confinement of a semiflexible chain is studied in the 
Odijk regime. It will be shown how as a consequence fluctuations in the confined 
directions get coupled. Also the influence of torque on the global persistence 
length is studied. An important result is a dramatic decrease of the elongation 
of DNA in nanochannels under moderate torques. 

The paper is organized as follows: In section 2 we set up the model and 
its Hamiltonian. We first consider an isotropic channel and show how the 
parametrization can be changed from a potential strength to a deflection length. 
We next proceed to treat the general anisotropic case. In section 3 we discuss 
some refinements, like the addition of a stretch- and twist-stretch-moduli, how 
to map the model to hard-walled channels and we briefly discuss how torque 
influences the elongation problem in nanochannels. 

2 The model 

The free energy of a freely rotating semiflexible chain in a confining harmonic 
potential was calculated in Ref. [2]. These calculations we now perform for a 
persistent chain under torsion. In fact confinement by a constant harmonic po- 
tential is usually not what is needed. In case of a hard-walled nanochanncl, the 
width of the channel is better approximated by a fixed standard deviation than 
a harmonic potential. The confining potentials in vivo and in single molecule 
experiments are usually also far from harmonic, while the standard deviation of 
the channel is used as a variational parameter [20, 3]. Nonetheless the harmonic 
potential is used as a Lagrange multiplier, an artificial confining potential to set 
the standard deviation <j of the transversal polymer distribution as measure of 
width of the confinement. The benefit is an analytically solvable model and as 
long as the confinement is strong enough specification of the second moment of 
the chain distribution will do. 

The typical length scale, that determines the physical properties is the 
deflection-length A = \jPb<J 2 of the chain, the length above which confinement 
dominates thermal fluctuations. The main calculations will be performed in 
a fixed number of turns, or linking number Lk, ensemble. To keep the treat- 
ment transparent we model the persistent chain as a framed space-curve with a 
persistence- length P& and a torsional persistence length P c . The results are eas- 
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ily extended with a finite stretch modulus S and a twist-stretch coupling, B [13]. 
All energies are scaled by ksT (unless mentioned) for notational reasons. As a 
result forces have the dimension of an inverse length. 

The chain is confined in two perpendicular directions, say x and y, by a 
harmonic potential. As in Ref. [2], we assume the confinement to be strong 
enough for the z-coordinate to be a single valued function of the curve. The 
framing is needed to describe the twist degree of freedom. Since the confinement 
potential is localized in space we cannot resort to the usual tangential description 
but have to start from the space curve [2]. The Hamiltonian is: 



t>L c 

J$?= I ds 



P b fd 2 r(s)\ 2 P c]2 



+ |rj(.) + |r»" 

As parametrization we choose the arc-length of the chain. The strength of the 
harmonic potential is set by b x ^ y . The twist angle ip( s ) is defined as follows: let 
t(,s) = r(,s) be the tangent and n(s) be one of the unit normals that define the 
local frame. The twist density is defined as the differential number of rotations 
of n around the tangential direction: tw(s) = j-(t A n) ■ n, which is a total 
differential and so we can define ip{s) := 27rtw(s), 

The twist density is not a small parameter and it is not independent of 
the radial fluctuations. To perform a perturbation expansion we split of the 
independent fluctuation degrees of freedom. This we can do using White's 
relation [21, 4], relating linking number, twist Tw, and writhe Wr (the Gauss 
linking number of a closed curve with itself): 

Tw := / ds tw = Lk - Wr (2) 
Jo 

By assumption the chain is almost fully elongated permitting us to take periodic 
boundary conditions. With this choice are the tangents at the ends of the chain 
parallel. Since 2tt times the writhe+1 of a closed curve is equal to the surface 
on the direction sphere enclosed by the path of the tangent mod 4-7T [5], and 
the tangents at both ends are parallel, we can define the writhe of the chain 
as the area enclosed by the tangent along the contour, implicitly setting the 
of the writhe to be the z-axis. Since the z-coordinate is single valued we can 
use Fuller's equation [5, 1] turning the writhe of the chain in an integral over a 
writhe density uj that is up to quadratic order in the fluctuations: 

wr ( s ) - -^{ f x{s)r y {s) - r x (s)f y (s)) (3) 

The linking number density we define as lk := Lk / L c . Now define a fluctuation 
field <f>(s) :=ip-2Trlk-2nwr(s). FromEq. (2) it follows that J^ c ds<f>(s) = 0. We 
use this to eliminate the twist term and can now integrate out the independent 
field </>. Making use of Eq. (3) results up to second order and an irrelevant 
constant in: 



l>L c 

J$?= I ds 



— r 2 4- —r 2 -A- hir 2 

2 2 x 2 y (4) 
irP c lk(f x r y - r x fy)} + 2ir 2 PJk 2 L c 
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The first term contains up to quadratic order only the x and y components since 
the z component of the tangent vector t(s) is quadratic in the fluctuations: 

r 2 (s) = i x {s) 2 + i y {s) 2 + i z (s) 2 ~ i x {sf + i y (s) 2 (5) 

In the following we will take r to be two-dimensional. For reasons of clarity we 
present the analysis in more detail for an isotropic confinement 



2.1 Isotropic confinement 

For isotropic confinement we have b := b x = b y . We can get rid of some of the 
coupling constants by appropriately scaling the contour lengths and the channel 
width [2] 

S = 6 1 /4p"V4 s f . = p fc l/8 6 3/8 r; (6) 

and so P c = 6 1 /4p" 1 / 4 p c but Ik = b^^P^lk. This scaling results in: 



,W = 



L, 





2<K 2 P r lkL r 



if 2 (s) + if 2 (s) - ^C(r x r y - f x r y ) 



with 



2TlP c lk 

P b 3/4 6i/4 



(7) 



(8) 



The configurational partition sum is a path integral over the two-dimensional 
field f : 

Z = J £>(r)e~^ (9) 

The last term of the Hamiltonian is constant and can be taken out from under 
the path integral. It represents the unperturbed twist energy. What is left is 
the fluctuation part which can be diagonalized, after a Fourier mode expansion 
using periodic boundary conditions, with eigenvalues: 

ev n ,± = pi + 1 ± Qpl with p n = ^ (10) 

L c 

Negative eigenvalues appear when £ > ( cr = 4/3 3 / 4 and so one can speak of 
a critical linking number density: lk cr = 26 1 / 4 P fc 3 ^ 4 /(3 3 / 4 7rP c ). The partition 
sum Z = J[Dr] cxp[~J^}, becomes a product of Gaussian integrals, that can be 
compared to the free chain [11]. The result can be written as: 

OO g 

Z = A/"cxp(-27r 2 p/fc 2 L c ) TT ^ (11) 

ev n + ev n ,- 

with Af collecting constant terms. In the denominator of the product we have 
a quartic expression in n 2 : Let r]i,fji(i = 1,2) be the pairs of conjugate roots: 

fa? + I) 2 - C 2 vf = 0. (12) 
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Then we can write the product as: 

oo oo 1 

II II II . 



i=l,2 n=l 



S~2 s ' m ( n V^) sin(7T^I) 



(13) 



The complex part of the square root of the rji's give an exponential decay with 
increasing chain length. Now write ^/?/i,2 = ^(#1,2 ± Wi,?)- In the long chain 
limit, as long as the yi are non-zero we obtain the free energy density as: 



f= lim -— logZ = 2tt a PJk 2 

L c — »-oc j_j r 



1/4 



Ph 



\Vi\ 



(14) 



The first term is the zero-temperature twist energy density stored in the chain, 
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Figure 1: Flow of the roots ^/rji in the complex plane under increasing torsion. 
The size of the arrows indicates the speed with which the roots flow as explained 
in the text. 

the second fth the thermal fluctuation contribution in presence of the confining 
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potential. In case the linking number density is zero, the roots are up to a sign 
equal: |j/i,2| = 1/V2 and we retrieve Burkhardt's result [2]. The flow of the 
roots under increasing linking number density is depicted in Fig. 1. Starting 
from Ik = the root pairs move apart. Both decrease their imaginary parts, 
decreasing the free energy. This does not happen in a symmetric way, as is 
shown by the line segments along trajectories, where each segment corresponds 
to a fixed increase of Ik, the arrow pointing in the direction of the flow. At 
lk cr the right pair of roots becomes real indicating a singularity in the partition 
sum. 

As explained above the harmonic potential is only there to set the confine- 
ment. The target is to express all quantities in the standard deviations of the 
transversal distributions of the chain in the channel. To do this and to extend 
the calculations to a non-isotropic confinement we perform the previous calcu- 
lations in an alternative way, taking first the logarithm of Eq. (11) and then the 
infinite chain limit: 



fth 



1 oc 



,V4 



Pb 

2nPJk 



(pt + I) 2 ~ C 2 Pl 
m dp 



o 



' log(£c 



2" 



log 



P 



(p 4 



AC 



V2- 



dp 
2^ 



log 1 



1)2 _ £2 p 6 

cv 



(15) 



(p 4 + l) 2 , 

expanding the logarithm, making use of the fact that (p 3 < p 4 + 1, for £ < £ cr , 
and noting that the series converges uniformly in any closed interval contained 
in [0,Ccr) this results in: 



2irPJk 



in 



p b c 

2ltP r lk 



71=1 J ° 



dp 
2^ 



2^6 



CV 



PbC 



with: 



(P 4 



n=0 



4n(2n)\ 



(16) 



This can be written as a linear combination of generalized hypergeometric func- 
tions: 



fth 



2-KPAk 



PbC 

3C 2 
+ T2 



V2 1 - 3 F 2 (-j, ^; j, |; ^r) 



—a p (1 _7_ 11. 3 3. C 

qo 3 2 V 4' 12' 12' 4' 2' Cf^ 



making the bifurcation point £ cr explicit. Defining: 



c 

x:= — 

S>cr 

results in a free energy density given by: 



:=V2 



2-KP c lk 



PbCc 



— - V a x 2 ' 1 " 1 = -- — - — + 0(x 7 



(17) 



(18) 



G 



From the series expansion of fth it is straightforward to find a series expansion 
for the standard deviation a r = J (r 2 ) /2 of the transversal position of the chain 
in the channel: 



dfth 
db 



Pbl 



OO 

n=0 



(2n — l)a n x 



2n+3 



(19) 



which is again a linear combination of generalized hypergeometric functions. Wc 
now apply Lagrange inversion to this scries in order to eliminate b in favor of oy 
as dependent parameter for the free energy density. To be able to use Lagrange 
inversion the derivative at the point of inversion should be non-zero, which 
requires non-zero coefficients of the linear term. The lowest term in Eq. (19) is 
cubic, so if we take the cubic root of this series we obtain a series to which one 
can apply Lagrange inversion. Introducing the deflection length, A := (P^of) 1 / 3 
we start from its series expansion: 



X(x) 



1 



/ao \ - r(- 

— x E- 



k) 



and apply Lagrange inversion to this: 

d n-l 



(20) 



x(X) 



E 1 

n=l 



! dx™" 1 \X(x) 



x=0 



The first terms of the resulting scries arc: 

7 3 A 3 7 7 5 A 5 



x(X) = 7 A 



6V3 



216 



G(( 7 A) 7 



f th (x) = l-^ + ^ + o iz y*) 

A 3^ 243\/3 



(21) 



or in terms of the unsealed variables 

1 



th 



1 - - ( nPJk— 
2 V P b 



24 I P b 



O 



nPJk— 
Pb 



(22) 



This free energy density contains the enthalpic free energy of the harmonic 
potential. But the harmonic potential is just a Lagrange multiplier to set the 
standard deviation, thus this enthalpic part of the confining potential should 
be subtracted in order to obtain the net cost for confining the molecule. The 
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cnthalpic contribution is easily calculated from the inverted series: 

fpot 



V 2 > = - 

2 W 16 



Pb 
X 



irPJk- 



Ph 



^(A) 
3 



TVPrlk- 



X 

Ph 



+0(iixP c lk— 



(23) 



Subtracting this free energy density from Eq. (22) gives the confinement free 
energy density: 



fconf — , 



nPJk^- 



nPJk^r 



8 



Ph 



Oi irPJk 



(24) 



Note that the free energy of confinement seems to decrease with increasing 
linking number, while the standard deviation is kept constant. This is actually 
the twist energy reduction due to thermal writhe. The unperturbed twist energy 
density was the last term in Eq. 7 adding the perturbation contribution we 
recover the effective twist free energy density: 

P r X 



ft 



2lT /l Pr 1 



Ik 



i ._ 



2^PJk z 



(25) 



where a "renormalized" torsional persistence length P c , is defined that is of the 
same form as that of a WLC under tension [15]. From Eq. (4) we see that 
the linking number density is conjugate to the thermal writhe density in f t h, 
allowing us to calculate its expectation value and from White's equation the 
expectation value of the twist density: 

. . 1 df th (lk,b) 3P c Xlk 



4n 2 P n 



(tw) ~ Ik — (wr) 



dlk 



3P C X 
W 



Ik 



(26) 
(27) 



Thermal motion decreases the extension of the chain. It is calculated by adding 
a tension term in the channel direction to the Hamiltonian: 



Fz~Jf+ I C dsF(l- \v 2 {s)). 
Jo 2 



(28) 



The extension in the z-direction is p := Az/L c = dpf(F)\F=Q- After a mode 
expansion the fluctuation part can be diagonalized again with F dependent 
eigenvalues. The calculation of the free energy density goes as before: 



fth 




(29) 
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For the contraction factor we obtain in terms of the deflection length: 



A 
2P b 



(30) 



A last quantity that is easily extracted is the expectation value of the torque r 
one has to apply to keep the linking number fixed: 

« = 2^tf = 2 * ft,t < 31 > 

The influence of torque on the confined polymer is small for a strong con- 
finement. The restricted space around the polymer hinders the production of 
thermal writhe. Upon increasing the standard deviation thermal withe becomes 
stronger. In the torsionless case, the demand that the z-coordinate is a single 
valued function along the backbone results in the demand that the deflection 
length is much smaller than the persistence length. With torsion, this upper 
bound decreases with increasing torque. From the series (29), (30) it can be 
seen that the demand on deflection length changes to 

Using Eq. (31) this can be written as a torque dependent constraint on the 
standard deviation: 

8/2 



"^U^^j < 33 > 

Note that r is scaled by fceT. In section 3 the effect of this approaching insta- 
bility on the extension is discussed. In Figure 2 the effect of a finite torque is 
shown on the semi-local extension of DNA, with torques that are realistic for 
torque producing molecular motors in vivo 1 . It is important to realize that the 
strength of the harmonic potential needed to obtain a fixed standard deviation 
increases with increasing torque. 



2.2 Non-isotropic confinement 

The strategy for the non-isotropic case will be the same as for the isotropic case. 
The increased complexity is a consequence of the twist mixing fluctuations in 
the two transversal directions. The Hamiltonian is again block diagonalizcd 
after a Fourier expansion. With 

= 33(6g6, + b x b\) + (g + 146.6, + b\fli - feg _ b 3 

1286.6, ' ^ ' 

and a scaling like in the isotropic case we see that the expression for the critical 
point is, like in the isotropic case: C CT = 4/3 3 / 4 . We assume £ to be small. 
Defining b xy := b x , y /b the fluctuation part of the free energy density can again 

1 see section 3 
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Figure 2: The local relative extension for DNA in a confined environment, 
producing a standard deviation up to 8 nm with a range of torques applied 
along the backbone. The persistence lengths were chosen to be Pi, = 50nm and 
P c = lOOnm 



be written as an integral over the determinants: 
\V* r ,oo 



fth = - 



Ph 



dp 
2^ 



log 



P 



™(*»_ r°° dp 

n L 9.7T 



27T \(p4 +bx )(p4 +by ) 



(35) 



Assume b x < b y . Integration results in a power series in £, with coefficients that 
can be expressed in series over hypergeometric functions. We will only need the 
first terms. They are: 



Jth 



1/4 (IT 



El/4 



V2 
rl/4rl/4 



+ b 



1/2 



(36) 



C 2 + 0(c 4 ) 



Since the terms in the series expansion are independent of the b— scale, we 
can choose it as b = l/P b 3 and change from b x>v to the dimensionless z X:V := 

(Pbb x , y ) 1/4 resulting in 



fth 



z x + z y z x z y( z x z x z y + z y ) 



Pb \ V2z x 



ty AV2( Zx + Z y)(4+ z2 y ) <:i + 



(37) 
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with Ci = 2irPJk. To get an expression for the free energy as function of the 
deflection lengths we will use a multivariate extension of Lagrange inversion [8, 
7]. As in the isotropic case we extract the standard deviations for the two 
channel directions from the free energy density: 

p3 

°l = — %-4 d zJth 

_ P b z l (-1 | z l z y(3 z l + 2z x z y + Zy) ^ 2 \ 
~ 2V2 \ 1+ ±(z x + z y y{zl + zlY C i+"- y )' 

and z x O z v for a y . With f x>y := z x _ y / X XtV (zx, zy), we have a relation between 
A^y and z XiJ/ of the right type to apply Good-Lagrangc inversion, that is: 

z x , y = fx, y K,y and f Xt y(0,0) ^ 0. (39) 

Under these condition given any (formal) power series g in z x and z y we can 
express g uniquely as a power series in X x , X y by the Good-Lagrange formula [ , 

7]: 

[X x . y ] n ' m g(z x (u x ), z y (u y )) 

= [Z X , y ] ' g( Z X, Z y)f x (Z X , Zy)fy (Z X ,Zy) 

Det(V * 

= [^, W ] n ' W 'a(^ ) Z y )f^(Z X ,Zy)fy n (Z X , Zy) 



L7 z. 9aj(z x ,z y ) 



^Det ~' . J ;-' tf/ , (40) 
9 V^^jZj) dzi / 

where [a;i j 2] n ' m -F(xi.2) denotes the coefficient of the monomial x\x™ in the 
(formal) power series F. When applied to fth we obtain a series that is symmetric 
against interchange of X x and X y . We define the symmetrized monomial: 



n,m]:=W + <W (41) 



For f tfl we get: 



f tt = - [-1,0] -^2 ([1,0] -[4,-3] 



+ [5,-4]- [8,-7] + ---)+o(c i 6 ) 
Ir_i nl # [3,i] + [2,2]/2 

2 l ' J 6P b 2 [1,0] [2,0] 

1 fA _i A _u _ 2(7rP c ^fc) 2 A%A„ + XjXl + A a Ag 

2 l * v ' 3P fc 2 ( A, + A y )(A2+A 2 ) 



(42) 



The combination of powers in the deflection lengths in the 2nd term will appear 
several times. For this reason we define an effective twist deflection length 

- X^tX v -(- A~Ar. -(- A^A?. 

A •= 9 f43^ 

' (X X +Xy)(Xl+Xl) ■ [ > 
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The factor 2 is a matter of convention, chosen for comparison with the P c 
rcnormalization concept of a chain under tension of Ref. [15]. Let X x be the 
smaller of the two deflection lengths, corresponding to the stronger confined 
direction, this effective deflection length has a value in between ^X x for \ x = X y 
and 2[1 - (X x /P b ) 3 ]X x for A^ << X y . 

Apparently the thermal writhe effect on the free energy is dominated by the 
smallest deflection length. This is understandable since writhe is non-planar and 
since the polymer is directed along one axes, tightening one of the transversal 
directions forces the polymer to be more planar. Just like in the isotropic case 
we subtract the free energy contribution of the harmonic potential, to get the 
pure confinement free energy. The calculations are merely a repetition of the 
fth calculation so we just give the result up to the same order: 



hot — o\ A x + A y ) H ^p2 A 



fconf = f th hot = + ^S^- ( 44 ) 

The second term can again be interpreted as a reduction of the twist free energy 
density. Adding it to the unperturbed twist free energy density we find: 

f twist = 2ir 2 P c (l - |j^J Ik 2 =: 2TT 2 PJk 2 , (45) 

the "renormalized" torsional persistence length P c , is again of the same form 
as the thermal softening of a WLC under tension [15]. The average writhe and 
twist densities follow from the same derivation as in the isotropic case, and using 
the Good-Lagrange formula for parameter change results in: 

, , XPJk , , / XP, 



AP 2 " ' V 4P h 2 



(tw) = 1 - — | Wfe (46) 



The extension is calculated like before by adding a force term (28). We perform 
the derivation by F before the integration, choose b = P 3 and expand in £i. 
The resulting extension is expressed in z 



4a/2 16V2 {z x + z v ){z 2 x + ziy 1 x ' 



Finally resorting again to Good-Lagrange inversion (40) we obtain 

1 / Aii- Ar 



P = 1 

4 V P b P b 



(7X 3 x X y -2X$X± + 7X x Xl)(<KP c lky 



(48) 



12(X x + X y )(Xl + X 2 y ) 2 p3 ■ 

Also this last expression has the feature that it is the tightest direction that 
determines the thermal writhe. In Figure 3 the extension of DNA upon non- 
isotropic confinement is shown. The amount of confinement anisotropy is changed 
by varying the standard deviations such that the sum of the two deflection 
lengths remains constant. In that way the contraction due to bending degrees 
of freedom does not change. The isotropic point is at a x = a y = 20nm. In 
the limit of vanishing standard deviation in one direction, thermal writhe is 
completely suppressed. 
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z = 1 pNnm 
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Figure 3: The squeezing out of thermal writhe for DNA by anisotropy for 3 
different torques while keeping the influence of bending modes contant. The 
elastic moduli are chosen as in Fig. 2 



3 Some extensions 

We did not include a twist-stretch coupling, that is typically nonzero for chi- 
ral molecules, or or a stretch modulus in the calculations. To include them 
complicates the calculations only modestly. Mainly because the stretch modu- 
lus is large compared to tensions normally applied, their influence is small. In 
lowest order it amounts to a correction of the bare torsional persistence length 
to P' c := P c — B 2 /S 1 where B is the (dimensionless) twist-stretch coupling, 
and S the stretch modulus [13]. Furthermore the contraction factor acquires a 
contribution: p — > p — 2tt Blk/S. 

The calculations were performed using a harmonic potential as a substitute. 
To relate this to a hardwalled channel we can extend the reasoning from Ref. [20] 
and calculate the deflection length of the hard-walled channel, with channel 
width d, by comparing the confinement free energy at zero linking number den- 
sity with the numerical results from Ref. [2, 22]. This leads to: 

One can reason now that since we have the free energy as function of the stan- 
dard deviations, and thus up to this numerical factor of the channel dimensions, 
we can extend this to a finite torque. Had we used the potential strength b as 
parameter the chain would fluctuate through the channel wall. 

A WLC under tension undergoes upon increasing linking number a buckling 
transition to a plectoneme, a critical multi-plectoneme phase or to a multitude 
of irregular configurations of loop- like structures (which might lead to structural 
changes in the WLC) [3]. Into which of these 3 categories the transition belongs 
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depends on material and environmental parameters, but the linking number 
density is in all cases controlled by a mechanical instability at a critical linking 
number density of lk cr = \J F / (?: P c ) . Comparing this with the instability in 
the confined chain, we see that the deflection lengths also here have a comparable 
role. The critical linking number density of the confined chain corresponds to 
that of a chain under tension F when: 



As an important example we mention DNA where the maximum linking number 
density causing a buckling transition is for forces of around 4 pN (F = 1 imr 1 ), 
higher forces causing a structural transition within the DNA double helix. With 
a persistence length of 50 nm this would correspond to a channel dimension 
~ 3 nm, which is an order of magnitude smaller than the channels used in 
confinement experiments. In a plcctoneme on the other hand the fluctuations 
of the strands can easily have comparable standard deviations [3] , showing that 
plectonemes are indeed stable. 



Figure 4: Extension of a 60 /im DNA chain in square nanochannels with widths 
from 30-450 nm, without torque and with torques of 2 and 4 pNnm (r = 0.5, 1), 
the dots are experimental data from Ref. [19] 

With increasing channel size hairpin formation starts to influence chain ex- 
tension [17] on top of local fluctuations. The free energy of an isolated hairpin 
in the mechanical limit, F mc is in the torsionless case that of an entropically 
squeezed mechanical hairpin [17]. The hairpins seen as rare defects in the lin- 
ear chain have a statistical density that increases exponentially with decreasing 
free energy of the hairpin. The global existence length is the average distance 
between hairpins When a torque is applied to a confined chain approaching the 
bifurcation point a local minimum will appear starting from an almost closed 









100 200 300 

Channel width in nm 



400 



14 



loop having a writhe close to one. In case the confinement is replaced by a 
tension, these loops form the nucleation point of plectonemes, super coiled he- 
lices perpendicular to the channel direction. The confinement does not allow 
for plectoneme formation and so the local minimum is the almost closed loop. 
These loops can not grow perpendicular to the channel axis, but can separate 
in two hairpins along the axes as shown in Fig. 5. It can be shown [3] that 




Figure 5: Two hairpins forming a planar loop of writhe ~ 1 

when the tangent of the ground-state varies on a much larger length scale than 
the thermal fluctuations, the writhe of the ground-state and thermal writhe 
can be treated separately, the total writhe being their sum. Conceptually the 
thermal writhe can be calculated from Fuller's equation with as reference curve 
the writhing ground-state. As long as this writhing ground-state is straight on 
the thermal fluctuations length scale the writhes separate. Since the confine- 
ment free energy (44) is only in higher order dependent on the linking number 
density, torque hardly influences the balance between entropic confinement and 
mechanical bending, that determines the hairpin free energy. The writhe of a 
pair of successive hairpins on the other hand does decrease the free energy and 
thereby the extension dramatically. In fact one can attribute a writhe of one 
half to a hairpin since two successive hairpins have this configuration as only 
minimum in their relative orientation, at least for a circular channel. Some 
minor contributions of other minima, reflecting the symmetry of the channel 
cross-section, might be present, but we will neglect them for the rest of the pa- 
per. Assuming the channel to be narrow enough that we can treat the hairpins 
as a dilute gas of defects we can follow the reasoning of Rcf. [17]. The free 
energy of a hairpin in the mechanical limit decreases by the work done by the 
torque: F mc (r) ~ F mc (0) — ttt. Or in terms of the global persistence length: 

g(r) = g(0)e-" T (50) 

As example Figure 4 shows the effect of moderate 2 pNnm (t = 0.5) and 4 pNnm 
(r = 1) torques on the elongation of a 60 ^m DNA chain in a square channel 
with widths ranging from 40 — 100 /im, following the theoretical treatment of 
Ref. [17] and in comparison with experiments from Reisner et al. [19] on DNA in 
nanochannels in the absence of torsion. The shortening by thermal fluctuations 
was calculated using equations (48), (31) and (49). As material parameters were 
used Pb = 57 nm (DNA with dye as in Ref. [19]) and P c = 100 nm. Applying a 
moderate torque induces its collapse mainly due to the lowering of the hairpin 
free energy. Only above 350 nm does the twist instability of confinement become 
important. A moderate torque of 2 pNnm reduces the elongation to less than 
half the value it had in the torqueless case for channels larger than 90 nm. 
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4 Summary and conclusion 



As is shown in this paper, torsion can have a large impact on the confinement 
of semiflexiblc polymers. For nanosizc confinement the thermal writhe caused 
by an applied torsion is damped by the tightest confinement direction. The 
apparent softening of the chain is similar to the softening of a semiflexiblc poly- 
mer under torsion [15] with the appearance of a new deflection length. The 
instability that appears in the free energy landscape sets new limits to the mag- 
nitude of the torque the molecule can propagate depending on the confinement 
of its environment. It is important that the strength of the confinement in one 
direction already suppresses the thermal writhe and thus increases the torque 
that propagates. 

An important tool in the investigation of biopolymers and their interaction 
with proteins is the possibility to spread out the molecule in nanochannels of 
diameter [19, 10, 23, 18] comparable to the molecules persistence length. The 
theory presented here shows that torques that for example RNA-polymcrase can 
induce [9] (up to 5 pN nm) do increase the demands on channel size considerably. 
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